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Abstract — We observe that in a multi-channel wireless system, 
an opportunistic channel/spectrum access scheme that solely 
focuses on channel quality sensing measured by received SNR 
may induce users to use channels that, while providing better 
signals, are more congested. Ultimately the notion of channel 
quality should include both the signal quality and the level of 
congestion, and a good multi-channel access scheme should take 
both into account in deciding which channel to use and when. 
Motivated by this, we focus on the congestion aspect and examine 
what type of dynamic channel switching schemes may result in 
the best system throughput performance. Specifically we derive 
the stability region of a multi-user multi-channel WLAN system 
and determine the throughput optimal channel switching scheme 
within a certain class of schemes. 

I. Introduction 

Advances in software defined radio in recent years have 
motivated numerous studies on building agile, channel-aware 
transceivers that are capable of sensing instantaneous channel 
quality 0, 0, 0. With this opportunity comes the challenge 
of making effective opportunistic channel access and trans- 
mission scheduling decisions, as well as designing supporting 
system architectures. There have been extensive studies on 
dynamic channel access in a multi-user, multi-channel wireless 
system, see e.g., 0, 0- By allowing users to dynamically 
select which channel to use for transmission, these schemes 
aim to improve the system performance, typically measured 
by the total (or per user) throughput, the average packet 
delay and etc, compared to a system with a single channel or 
more static channel allocations. The main reason behind such 
improvement lies in temporal, spatial and spectral diversity. 
That is, the quality of a channel perceived by a user is time- 
varying, user-dependent, and channel-dependent. 

Within this context we make the additional observation that 
there is also a congestion diversity in that a channel with 
fewer number of competing users presents better quality for 
a user. This is particularly true in a random access setting, 
where a large number of competing users can induce large 
backoff timer values that in turn lead to longer waiting time 
and lower throughput. We note that in a multi-channel system, 
an opportunistic access scheme that solely focuses on channel 
quality sensing as a result of random fading and shadowing, 
e.g., by measuring received SNR 0, 0, may induce users 
to use channels that, while providing better signals, are more 
congested. This can reduce the expected performance gain, 
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or even turn gain to loss. Ultimately the notion of "channel 
quality" should include both the signal quality and the level of 
congestion, and a good multi-channel access scheme should 
take both into account in deciding which channel to use and 
when. 

Motivated by the above, in this study we examine the 
possibility of utilizing congestion diversity to promote certain 
performance measures, e.g., throughput. As mentioned above, 
our ultimate goal is to construct an opportunistic channel 
access scheme that is aware of both signal quality and user 
congestion. However, in the present paper we will limit our 
attention to addressing the congestion aspect only; a good 
understanding of this aspect is a crucial first step in this effort. 

Specifically, we ask the question of what type of dynamic 
channel switching schemes will give the best performance in a 
multi-channel WLAN. This will be evaluated using the notion 
of stability region of a scheme. This is because more effective 
resource allocation and sharing can achieve a lower overall 
congestion level, thus expanding the range of sustainable 
arrival rates and resulting in a larger stability region. The 
scheme with the largest such region is commonly known as 
the throughput optimal scheme. With this objective, we set 
out to study the stability region of a multi-channel WLAN 
system where users are allowed to dynamically switch between 
channels. 

The main contributions of this paper are as follows. 

• A mean-field based model is constructed to characterize 
the stability region of a multi-channel WLAN system. We 
show that the size of the backoff window plays a decisive 
role in shaping the corresponding stability region: when 
the backoff window is sufficiently large, the stability 
region is convex; as the window size decreases it evolves 
into a concave region. 

• Using this mean-field model, we provide an analyti- 
cal justification of using channel- switching policies that 
achieve load balance in systems with symmetric chan- 
nels. This is then extended to systems with asymmetric 
channels. 

• We propose several simple heuristic implementations of 
the channel- switching policies presented in this paper. 

802.11 DCF has been very extensively studied in the liter- 
ature, ranging from throughput performance in the saturated 
regime 0, and the non-saturated regime 0, fTOlh to its 
rate region ifTTIl . fT2ll . to name a few. To the best of our 
knowledge, however, none has studied multi-channel WLAN 
in the context of stability region. Works most relevant to ours 



include ones on the stability region of slotted Aloha (e.g., f\3l ) 
and the rate region of 802.11 DCF fTfl, P. 

In the remainder of the paper, we first introduce a system 
of equations to characterize the stability region of a single 
channel WLAN consisting of multiple users within a single 
interference domain (SectionllTfl) followed by numerical results 
(Section HVl). We then extend the same method to characterize 
the stability region of a multi-channel system and use this 
result to determine the throughput optimal channel switching 
schemes within a class (Section [VT). We also discuss how such 
schemes may be implemented in practice (Section IVI1). 

II. System Model and Preliminaries 

Consider a multiple access system using the IEEE 802.11 
DCF. There are TV nodes (or users interchangeably), indexed 
by the set N = {1,2, ...,iV}, each with an infinite buffer, one 
transceiver and uses the same parameterization. We assume 
the channel is ideal and there is no MAC-level packet discard, 
i.e., there is no retransmission limit of a packet after collision. 
Throughout the analysis we also adopt a few other simplifying 
assumptions to make the problem tractable; these will be stated 
in the context to which they apply. It should be noted that due 
to the complexity of the problem, successive simplification in 
the modeling effort is a rather common practice and has been 
used in most if not all previous works. We later show that 
these simplifications do not impact the accuracy of the model 
under normal operating parameter values. 

The key to our method is to model the queue at each node 
with a service process defined by 802.11 DCF as a slotted 
mean field Markov chain [ITU . 

Definition 1: A virtual backoff timer of the system (or of a 
virtual node) is a universal timer for all nodes in the system: it 
counts down indefinitely, alternating between the count-down 
mode (when nodes in the system are counting down) and the 
freezing mode (when some node in the system is transmitting). 
The slot time is thus a random variable. 

Remark 1: The above definition provides a universal slot 
time for all nodes in the system, and we will assume that the 
real backoff timer at each node is synchronized to this virtual 
timer on slot boundaries. The motivation behind such a con- 
struction originates from the principal difficulty in modeling a 
non-saturated system: the service process at each node runs in 
embedded time in terms of a slot, which is a random variable, 
whereas the packet arrival process is more naturally described 
in real time fT4ll . This difficulty does not exist in saturated 
analysis, see e.g., Q, where arrival processes do not play a 
role. 

We next introduce three key assumptions in our model, 
followed by a discussion on their implications and limitations. 

(Al) The MAC layer arrival process at node i is Poisson with 

rate A^ bits per second. 
(A2) The service time of a packet, i.e. the time from the initial 

backoff to successful transmission, is (i) exponential with 

service rate pi at node i, and (ii) independent of all arrival 

processes. 



(A3) Let S(t) be the counting process of the number of slots 
accumulated up to time t. S(t) is assumed to be (i) 
independent of Qi(t), and (ii) renewal. 

Denote by A = (Ai, A2, . . . , Xn), and by {Qi(t)}t the 
queueing process at node i (also written as Qi(t) for sim- 
plicity), i.e., the number of packets in node z's MAC queue at 
time t. We now formally define the stability region of system 
as follows. 

Definition 2: The stability region A is the set of all A G 
such that Qi(t) admits a stationary distribution for all i with 
arrival rates A under the 802.11 DCF scheme. 

The above simplifying assumptions are not entirely realistic. 
Typically, due to congestion control by upper-layer protocols, 
e.g., TCP, the arrival process to the MAC layer is neither 
Poisson nor independent of the service process. However, as 
our objective is to explore the inherent properties of 802.11 
DCF, the independence assumption is adopted to decouple 
the MAC layer from upper layers, while the Poisson and 
exponential assumptions are adopted to avoid technicalities 
that can obscure the main insight. Note that under the mean 
field methodology, each node is analyzed in isolation from the 
activities of all other nodes which are collectively regarded 
as an aggregate stationary process. Within such a framework 
the packet service time is taken to be stationary (see e.g., 
Bianchi's well-known mean field Markovian model of the 
service process Q). 

With Al and A2, Qi(t) is then a well-defined M/M/l 
queue, and for a given A, A G A if and only if Qi(t) is positive 
recurrent. Equivalently we may consider the utilization factor 
pi at node i, given by pi = minj^-, 1}: the queue is stable 
if and only if pi < 1. If Qi(t) is positive recurrent, then it is 
ergodic and we have lim^oo P(Qi(t) > 0) = 1 — 7^(0) = pi, 
where m is the stationary distribution of Qi(t). If Qi{t) is 
transient or null recurrent, in which case pi = 1, we have 
limt^oo P{Qi(t) = 0) = = 1 — pi. Therefore, pi is 
asymptotically given by lim^oo P(Qi(t) > 0) in all cases 
in our model. 

For technical reasons we will also consider the embedded 
queueing process Qi(n), n = 1, 2, • • • , defined as Qi(n) := 
Qi(T n ), where T n is the time of the nth slot boundary. Qi(n) 
is thus a discrete-time process constructed by observing Qi(t) 
at slot boundaries. 

For an arbitrary process S(t), Qi(n) is not necessarily 
Markovian. However, given assumption A3, durations between 
slot boundaries are i.i.d., constituting sampling periods that are 
independent of Qi(t). Hence Qi(t) is a discrete-time Markov 
chain under our assumption^. It's worth noting that A3 does 
not exactly hold in reality because the slot length is a function 
of a node's activity, and thus the state of its queue, even 
with the mean field simplification of other nodes' behavior 
(this is more precisely shown in the appendix). However, this 
dependence weakens when the number of nodes or the backoff 
window size is sufficiently large. We empirically show that this 
assumption does not impact the accuracy of prediction even 
with a small node population and backoff window size. 

1 This claim is proved in Appendix [G] 



Let pi denote the utilization factor under the discrete- 
time system Qi(n). In general pi ^ pi. Indeed we show in 
Appendix [A] that pi < pi where equality holds if and only if 
pi = 1 or pi = 0, i.e., node i is either saturated or idle. Similar 
to pi, pi is asymptotically given by lim^oo P(Qi(n) > 0). 

We will adopt Bianchi's decoupling approximation Q as 
another key assumption, stated as follows. Define Ci(j) := 1 if 
the jth attempt by node i results in a collision, and Ci(j) := 
if it results in a success. 

(A4) [Bianchi's Decoupling Approximation] For each node 
i G J\T, the collision sequence {Ci(j)} is i.i.d. with 
P(Ci(j) = 1) = pi for some constant pi. 

In reality successive attempts by the same node may occur 
if it repeatedly selects timer value while other nodes' timers 
remain frozen. In such cases the above assumption ceases to 
hold. This phenomenon can be prominent when the window 
size is small, and has been taken into account in some recent 
work fT5l . We call the string of successive attempts a run 
of attempts, and the first attempt in a run a run- first- attempt 
or simply first-attempt. If successive attempts were consider, 
we could alternatively adopt the following assumption: let 
C[(j) := 1 if the first- attempt of the jth run of attempts 
by node i results in a collision, and C[(j) := if it results 
in a success, and assume the first-attempt collision sequence 
(FACS) decoupling, i.e., 

(A4') [FACS Decoupling Approximation] For each node i G Af, 
the first-attempt collision sequence {C-(j)} is i.i.d. with 

P(C-(j) = 1) = Pi for some constant pi. 

In this study we will ignore the possibility of successive 
attempts for simplicity of presentation and adopt (A4). This 
is reasonable when the initial window size is sufficiently 
large. Our empirical results are fairly close between with and 
without consideration of successive attempts for large backoff 
windows. For small backoff windows, the discrepancy between 
the two will be illustrated in the numerical results. 

We will use the term backoff length to mean the total number 
of slots that a node spends between two successive timer 
renewals during the service process, which is the selected timer 
value plus 1. Define Nf and Nf x , respectively, as the numbers 
of slots and transmission attempts that node i takes in serving 
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We next derive a relationship between the transmission 
attempt probability and pi. Let r^(n) be the probability that 
node i initiates a transmission attempt in the nth slot. 

Lemma 1: n := lim^oo r^(n) exists and is given by Ti = 
k/Wi. 

Proof: Let Tx(n) denote the event that node i initiates 
an attempt in the nth slot. Then 

n (n) = P(Tx(n)\Qi(n) > 0) • P(ft(n) > 0) + 
+ P(T x (n)\Q t (n) = 0) • P(Q % (n) = 0). 

Consider now the sequence of slots in which node i has a 
packet in service. Given the decoupling among nodes, the 
occurrences of slots in which node i starts the service for a 
packet thus form renewal events. Regarding each transmission 
attempt as one-unit reward and using the renewal reward 
theory, we then obtain 



lim P(T x (n)\Qi(n) > 0) 
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Since lim^oo P(Qi(n) > 0) = p u and P(T x (ri)\Qi(ri) = 
0) = 0, the result follows. ■ 

To put the above result in context, one easily verifies that in 
the extreme case where all nodes are saturated and identical, 
we have pi = pi = p = 1 and Pi = p for all i. Consequently, 
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which is exactly the same as obtained in [7] Eqn (7). 

III. Single Channel Stability Region 

A. The stability region equation E 

Our first main result is the following theorem on the quanti- 
tative description of A. Let E[5^ q denote the conditional 
average length of a slot given that the queue at node i is 
non-empty but i does not transmit in this slot. T s and T c 
denote the lengths of a successful transmission and a collision, 
respectively. 

Theorem 1: A G A if and only if there exists at least one 
solution r = (ti, T2, . . . , tn) to the following constrained 



where W is the size of the initial backoff window and m is 
the value of the maximum backoff stage. Also note E[A^ X ] = 



2 We note that 1161 used a similar technique in computing the conditional 
transmission probability defined therein. 
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system of equations (E, C, A): 

Pi 
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subject to 
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(c) 

(i) 
(ii) 



< pi < l,Vi 

where P is the packet pay load size. 

Proof: Etfa)) is the result of Lemma [T] and is an 
immediate consequence of the definition of pi . Let the average 
service time at node i be X{ seconds per bit; the average 
service time per packet is thus PXi. Define Yi(j) as 
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Physically, is the average time between the beginning 

of the jth transmission attempt, which results in a collision, 
and the the beginning of the (j + l)th attempt, given that node 
i encounters at least j collisions before completing the service 
of some packet. Since the collision sequence is geometric, we 
have 



W + l 



k=0 

X (pi) k (l- Pi ) 



oo oo 



3=1 k=j 

W + l 



W + l 



1 x 



J i,Q,T x l 1 

oo 

Therefore, 



1 E[5, 



i,Q,T x 



3=1 L 



2min{j,m}^ _ 



2min{j,m}^ _|_ j 



1 X 



= E 



E[S. 



1- Pi 1 *^* J ° 

Note that < 1 for all z, and we have ^ < 1 for all i 
as a result. In addition, K[S i q^J is finite (computed in the 
appendix). Hence we conclude that the packet service time 



is finite. Thus, the utilization factor of node i is given by 
pi = mmjA^X^, 1} and Etfc]) follows. C(0) is for the validity of 
r as a probability measure. (E, C(|i]), A) then constitutes a full 
description on the system utilization. 0© is the necessary and 
sufficient condition for stability as commented in the previous 
section. ■ 

For a given set of system parameter values, two sets of 
quantities are needed to compute E: E[5^ ojf J and pi,\/i 6 
A/\ These are computed in Appendix |B] and respectively. In 
particular, in AppendixOwe show that though it is analytically 
intractable, pi is well approximated by 



pi 



ftE[S.^] + (l-pi)E[Si,, 



where E[5^q] (resp. E[S^ q]) is the conditional average length 
of a slot given that the queue at node i is non-empty (resp. 
empty) at the beginning of this slot. 

B. Characterizing the solutions to E 

Without the stability constraint 0©, (E, C(Q), A) can be 
rewritten as a vector equation in [0,1]^, r = T(r), where 
t — ( r i? r 2, • • • , ^at) G [0, 1]^, and the existence of solutions 
can be shown by Brouwer's fixed point theorem. However, 
the uniqueness of its solution is in general difficult to prove; 
nevertheless, under the condition of a sufficiently large initial 
backoff window W, we have the following result on the 
uniqueness of its solution. 

With a large initial backoff window W, the probability of 
collision is small, so we have Wi ~ W ^ • We a ^ so observe 
that E[S»,q] « E[S^g] when W is large (cf. Appendix ©. 
Consequently, we can approximate by p^. Also, using the 
first-order Taylor approximation, we have Yij^i T^F ~ 1 + 
Sj/i r i f° r sma ll r - Let T s = T c = T for simplicity of 
presentation. Then E can be approximated by the following 
system of equations, 



Pi 
w+i 



(a) 



E: i 



Pi = 



W - 1 



T 1 



(b) 



Proposition 1: ((E), A) admits a unique solution. 

Proof: See Appendix El ■ 

Remark 2: 1) The above result suggests that E has a 
unique solution when W, the initial window size, is sufficient 
large. As an approximation we will take this condition to be 
equivalent to a large average backoff window. This is because 
the probability of a (first-attempt) collision decays inverse- 
linearly in W, and thus Wi is dominated by W when W is 
sufficiently large. 

2) As we will see numerically in the next section, multiple 
fixed point solutions may arise when W is small; this will be 
referred to as multi-equilibrium (as opposed to "multistable" 
or "metastable" iTUfll to avoid confusion). 

In the proof of Proposition Q] we in fact obtained the ap- 
proximated unique solution to (E, A). Therefore, by imposing 
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(c) W = 32, m = 5 
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(d) W = 128, m = 5 



Fig. 1. Solution components for various scenarios: an illustration. 



feasibility constraints C, we can induce a simplified version 
of (E, C, A) which is an approximation to A and is easier to 
compute. 

Corollary 1: When W is sufficiently large, A is approxi- 
mated by 

where 7^) = ±f / (l + and 7 ?(A<) = 

p(w+i) /v ^ 

Within the context of a unique solution to (E, C, A), con- 
sider A as input parameters and rewrite E as F(r, A) = 0, 
with (n + n) unknowns (r^'s and A^'s). We can then inspect 
the existence of an implicit function of r in terms of A, and for 
this we need to examine the invertibility of the corresponding 
Jacobian matrix. Note also that the correspondence between 
pi and (A, r) given by Eflc]) is a continuous function. If the 
Jacobian is invertible on the boundary of the stability region 
A in the space R+, then the continuity of pi = Pi(X) is 
established. Hence, on the boundary of A, denoted by dA, 
there exists at least one node i such that pi = 1. Unfortunately, 
the invertibility of the Jacobian on dA is highly non-trivial 
to determine and in general analytically intractable when the 
number of nodes is large. In the next section we numerically 
evaluate (E, C, A) and more is discussed. 

IV. Numerical Results: Single Channel 

Using (E, C, A), we can quantitatively describe the stability 
region of a single channel system, and some numerical results 
for the two-user case are illustrated in this section. The 
parameters used in both the numerical computation and the 
simulation are reported in Table U in Appendix [G] Under the 
basic access mechanism of DCF we have 

T s = + Header + ACK + DIFS + SIFS + 26 

where 6 is the propagation delay. 

A. Multi- equilibrium and discontinuity in p 

We first illustrate the existence of multi-equilibrium solu- 
tions and discontinuity of pi{\) in A; this is shown in Figure[T] 
We fix the value of A2 and increase Ai from to 4.5 Mbps. 



For each pair A = (Ai, A2), we solve for the fixed point(s) of 
E with the same set of initial values of n and pi for i G J\T 
to which we refer as a set of initial conditions (ICs). We then 
convert the results to p = (pi,i G A/") using Eqn. E©. The 
collection of the pairs (A, p(A)) then constitutes a solution 
component for this set of ICs. Note that this is obtained 
by solving (E, C(|i]), A) without considering the stability 
constraint 0©. We repeat the above computation for different 
sets of ICs under the same system parameters including W and 
m. The entire process is then repeated for different pairs (W, 
m). For each pair (W, m), the resulting solution components 
constitute an overall correspondence between the vectors A 
and p(A), and this is plotted for pi vs. Ai in Figure [T] 

In the first scenario as shown in Figure |l(a)[ where the 
initial window is of the smallest possible size for two users 
and window expansion is disallowed (m = 0), three different 
zones of the correspondence pi(Xi) are present, labeled as A, 
A' and B in the figure. In zones A and A\ a single fixed 
point is admitted and pi(Ai) reduces to a function, while in 
zone B we see two solutions. Along each solution component, 
there is a jump in pi in zone B as Ai increases; this is 
essentially a phase transition from stable to unstable regions. 
What this result illustrates is that depending on the initial 
condition, certain input rates may or may not lead to a feasible 
solution (a point in the stability region). Thus when such multi- 
equilibrium exists, we may have a collection of stability region 
A's given different initial conditions, and this phenomenon is 
illustrated in Figure [3] and discussed in the next subsection 
in detail. Recall that under our definition of stability region 
and Theorem [T] an arrival rate vector is considered within the 
stability region as long as there exists such an initial condition 
that induces so; the stability region thus defined is therefore 
the supremum of this collection when multiple equilibria exist. 
The advantage of the "stability region" A is that the points 
within are stable independent of the initial condition. With a 
slight abuse of terminology, we would later refer to this region 
as the stability region with multi-equilibrium. 

Intuitively, initial conditions with large values suggest a 
pessimistic prediction on the system stability under A, and 
it may thus result in a small A; by contrast, ICs with small 
values render an optimistic one and a larger A. Empirically, 
we find that the set of ICs with Ti = pi w 1 for i G M results 
in the earliest jump in pi and the one with n = pi = for 
i G N gives the latest. Consequently, solution components 



6 



Numerical boundary of A 

-b -Simulated boundary of A 




Numerical Boundary of A 

-B- Simulated Boundary of A 



1 2 3 4 5 

Ai (Mbps) 

(a) W = 128, m = 5 
Fig. 2. The stability regions in various scenarios - part I. 
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(a) Using Bianchi's decoupling approxima- (b) Using the first-attempt decoupling ap- 
tion proximation 



Fig. 3. The stability regions in various scenarios - part II: W = 2 and m = 0. 



resulting from these two sets of ICs define the boundary of 
zone B and the corresponding stability regions, forming the 
empirical supremum and infimum of the collection of A's. 

Inspecting the set of figures Fig. |l(a)||l(d)| we see that as 
the initial window increases, the multi-equilibrium gradually 
vanishes and the gap in pi caused by the jump discontinuity 
closes. 

B. Numerical and empirical stability regions 

We numerically solve (£, C, A) with two nodes to obtain 
the corresponding A, and then compare it with the simulated 
boundary. In simulation, for each fixed A2, we increase Ai 
with a step size A A, and compute the empirical throughput of 
node i obtained under A, denoted as S* 9 and the number of 
backlogged packets at node i by the end of simulation, denoted 
as B*. The simulator declares a point A unstable if there exists 
at least one i such that < A; and B^P/(XiTf) > ai, by the 
simulation time Tf, where is an instability threshold, < 
ai < 1. In the experiment we set A A = 0.1 Mbps (100 Kbps), 
Tf = 10 sec and = a = 1%. The stable point (Ai,A2) 
such that (Ai + AA,A2) is unstable is declared a point on 
the simulated boundary; the experiment is repeated for each 
A2 and the empirical mean value of Ai is recorded. Due to 
symmetry, only half of the boundary points are evaluated. The 
results are shown in Figure [2] 

Our main observation is that when the initial (or average) 
backoff window is large, the stability region is convex (Figure 
[2a)). The convexity gradually disappears as the window size 



decreases and the region is given by a near-linear boundary 
in Figure Etb). It becomes clearly concave when the window 
size is small (Figure Etc)). Interestingly, the case of W = 32 
is the most frequently studied in the literature, and a linear 
boundary of the capacity region has been observed in ifTTIl . 
As shown here, this linear boundary is only a special case in 
a spectrum of convex-concave boundaries. It is worth noting 
that in 021, Leith et al established the general log-convexity 
of the rate region of 802.11 WLANs. This implies that the 
rate region could be either convex or concave, though lTT2l did 
not associate this with the window size as we have explicitly 
done here. It also suggests that the rate region and the stability 
region may be quite similar in nature; this however is not a 
formally proven statement, nor are we aware of such in the 
case of 802.11. 

The change in the shape of the stability region as W changes 
may be explained as follows. Small W represents a highly 
aggressive configuration. This is much more beneficial when 
there is a high degree of asymmetry between the users' arrival 
rates. This is reflected in the concave shape of the region. 
When W is large, users are non-aggressive, which is more 
beneficial when arrival rates are similar, resulting in the convex 
shape. Numerically, the W = 8 case gives the largest stability 
region. This seems to suggest that the largest stability region is 
given by the smallest choice of W such that a unique feasible 
solution to (£, C, A) exists. It would be very interesting to see 
if this could be established rigorously. 

In Figure [3l we compute the stability regions of the case 
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where W = 2 and m = for two different sets of ICs. As 
discussed earlier, when multi-equilibrium exists we may have 
a collection of stability regions. This is clearly seen in Figure 
[3j three different zones A, A' and B in the correspondence 
pi(Ai) are mapped accordingly onto A. From these results, we 
may interpret that in zones A (A'), the system is uniformly 
stable (resp. unstable) regardless of the IC, while in zone 
B the stability of system depends on the IC. As noted in 
|[T4l . the simulated boundary reflects time-averages of multiple 
equilibria. 

As mentioned earlier, for small backoff windows the occur- 
rence of successive attempts is non-trivial, which our model 
has ignored. The first- attempt decoupling approximation A4' 
mentioned after A4 captures the nodal behavior more accu- 
rately, and the adaptation of E using this alternative assump- 
tion reduces to the computation of (conditional) mean slot 
length and is detailed in Appendix |B] In Figure [2b), we 
plot the counterpart of Figure Oa) using the first-attempt de- 
coupling approximation, and the discrepancy between results 
obtained using these two assumptions does exists. This is most 
notably shown in the numerical boundary A. The fact that 
the simulated boundary is now in between the two numerical 
boundaries verifies that this alternative assumption is more 
accurate. We do note however that for large windows this gap 
diminishes judging from numerical observation, which is to 
be expected. 

C. Discussion: from 802.11 DCF back to Aloha 

We next recall results on the stability region of slotted 
Aloha, the natural prototype of modern 802.11 DCF, and 
provide an intuitive argument on why the qualitative properties 
of the stability region shown in the previous section are to be 
expected. 

In ifTTlL Massey and Mathys studied an information the- 
oretical model of multiaccess channel which shares several 
fundamental features with slotted Aloha. They investigated the 
Shannon capacity region of this channel with n users, which 
is shown to be the following subset of \ 



C 



vect ( 71 ]^[(1 



0<Ti<l,l<i<7l 



where vect(^) = (vi, i>2, • • • , v n ), and ti is the transmission 
attempt rate of user i. In |[T3lL Anantharam showed that the 
closure of the stability region of slotted Aloha is also given by 
C, under a geometrically distributed aggregate arrival process 
with parameter 1 / A^) and probability A^ / . Xj that such 
an arrival is at node i. 

The above result on slotted Aloha can be used to explain the 
stability region of 802.11 DCF. Note that the main difference 
between the two lies in the collision avoidance mechanism. 
Instead of attempting transmission with probability < p < 1 
in a slot under slotted Aloha, under DCF each user randomly 
chooses a backoff timer value within a window. The effect the 
average backoff length W has on transmission under DCF is 
akin to that of restricting the attempt rate p within an upper 
bound = under slotted Aloha. Hence, the stability region of 
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Fig. 4. The stability region of slotted ALOHA and induced subsets. 



802.11 DCF may be viewed as a subset of C provided that 
we properly scale a slot to real time. 

To verify this intuition, let C w be the subset of C_when 
< Pi < i for all i. In Figure H we plot C and C w with 

different values of W. We see that as W grows, C w evolves 
from a concave set to a convex set, consistent with what we 
observed of 802.11 DCF in the previous subsection. It must 
be pointed out that this connection, while intuitive, is not a 
precise one technically. For instance, this connection might 
suggest that the stability region of 802.11 DCF will reduce to 
C when the average backoff length is 1 . This is however not 
true. In this trivial case, the stability region of 802.11 DCF 
is reduced to one dimension, i.e., the system is unstable for 
n > 2. This is because the retransmission probability of DCF 
is also lower bounded by the reciprocal of the window size at 
its backoff stage, and in the case when the backoff length is 
one another collision occurs with certainty. 

V. Multi-channel Analysis 

Using a similar, mean-field Markovian model as we did in 
the single channel case, we can show that the stability region 
of a multi-channel system under a certain switching policy g 
is given by another system of equations denoted as (E g , C, A), 
under the arrival rates A = (Ai, A2, . . . , Ajv), and subject to 
the feasibility constraints C; this is given later in the section. 
In addition to the same set of assumptions made in the single 
channel model, we assume that the system has K channels, 
indexed by the set C = {1, 2, . . . , K}. 

The fundamental conceptual issue accompanying channel- 
ization is the notion of a channel switching policy, either cen- 
tralized or distributed, that introduces channel occupancy and 
packet assignment distributions for each node. An additional 
technical issue induced by channelization is the heterogeneity 
of embedded time units among different channels. Since the 
slot length in a channel is by nature a random variable that 
depends on random packet arrivals, channels are in general 
strongly asynchronous in the embedded time units. Thus, as 
nodes switch among channels, we may need to switch the 
corresponding reference of embedded time in the slot based 
analysis. We therefore define the notion of a slot in different 
contexts as follows. 
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Definition 3: Consider the virtual backoff timer defined 
earlier separately for a single channel. A channel-slot (c- 
slot) is defined as the time interval between two consecutive 
decrements on this virtual timer for a given channel. 

Definition 4: Consider a virtual backoff timer at each node 
that counts down indefinitely according to the node's backoff 
state, and is synchronized to the virtual timer of the channel 
in which the node resides and is done upon switching. A 
node-slot (n-slot) is defined as the time interval between two 
consecutive decrements on a given node's virtual backoff 
timer. 

Remark 3: There is no inherent difference between the two 
types of slots. However, this differentiation of time refer- 
ences becomes crucial when we define quantities based on 
the random embedded time. This observation will be made 
more concrete in the analysis. We will also omit the explicit 
association of a channel (node) index with a slot whenever it 
does not cause ambiguity. 

A channel switching or scheduling policy g induces a 
number of distributions related to E g . Denote by Q™ (j) = 
{q\ k \j),k G C}, where q\ k \j) is the probability that node i 
is in channel k at the beginning of its jth n-slot, tj. Qf(j) is 
referred to as the the channel occupancy distribution in n-slots 
of node i in the jth n-slot. 

Denote by Q^j) = {qf\j),k G C}, where qf\j) is the 
probability that node i is in channel k at the beginning of its 
jth c-slot, ij. Qi(j) is referred to as the channel occupancy 
profile of node i at the jth c-slot. Note that Q£(j) is not 
necessarily a distribution and ^Z keC qt\j) need not be 1 for 
a given j. 

Denote by Qf(Z) = {qf\l),k G C}, where qf\l) is the 
probability that the Ith packet of node i is served in channel k, 
and Qf (Z) is referred to as the packet assignment distribution 
of node i. 

We have the following assumptions on policy g. 

(A5) Under g, Bianchi's approximation is still satisfied. 
(A6) g is independent of the binary state of the queue at any 

node (empty vs. non-empty). 
(A7) g is nonpreemptive in a channel for the entire service 

process of a packet; that is, a channel- switching decision 

is only made before or after the service process of a 

packet. 

(A8) The limits of Q^(j), Q-(j) and Q?(Z) exist under g 
as their respective arguments tend to infinity, and are 
denoted by Qf , Q\ and , respectively!! 

Similar as in single channel analysis, we impose the 
Markovian assumption on the discrete-time queueing process 
Q[ k \n), which is the embedded process of Qi(t) (queue state 
of node i) sampled at the boundaries of c- slots of channel k, 
and define pf ] = lim^^ P(Qf\n) > 0). Also, let rf } (n) 
be the probability that node i initiates a transmission attempt 
in the nth c-slot of channel k. Then we have the following 

3 These limiting quantities are related by well-define correspondences, 
which are detailed in Appendix [D] and those relations are used to numerically 
evaluate the stability region equation for a multi-channel system presented in 
this section. 



lemma; its proof is similar to that of Lemma \T\ (based on A6 
and A8) and omitted. 



Lemma 2: r 



:= lim n ^ oo r- k \n) exists and is given 

u (k) ^(k) Jk) iTffW u TT7( fc ) E[7V ? s ' (fc) ] . , u 

by r\ } = q\ >p\ ) /W i , where W i := E[ ^, (fc) J ] is the 
average backoff length of node i in channel k, with JV?'^ 

tx ( /c) 

and N i ,y ' defined in parallel as in the single channel case. 

(k) ^ 



Remark 4: Under A7, W i is given by 
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|(l- P f ) )^(2pfy 
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where p\ ' is the probability of collision in channel k given 
a transmission attempt and W is the initial backoff window 
size. 

Given any scheduling policy g, let A g be the corresponding 
stability region, and we have the following theorem character- 
izing A g . 

Theorem 2: A G A § if and only if there exists at least one 
solution r = (r W ,keC) where r W = (r^ k) , i G AT) to the 
following constrained system of equations (H g , C, A), 



Ak) 



4k) ,(k) 
Qi Pi 

r(k) 



(a) 
(b) 



keC 
(fe) 



w < -1 F r«? fe i 



■Pi 



T ( k) \ 



subject to 



C : 



< r\ k) < 1, Vi, k 
< pi < 1, Mi 



(c) 



(i) 
(ii) 



where i € Af and k € C; P is the packet payload size; 
E[5^^ ] is the conditional average length of a c-slot in 
channel k given that the queue at node i is non-empty but 
i does not transmit in this slot. 

Proof: The proof is an immediate extension of the proof 
of Theorem 1, given assumptions on g. ■ 

The existence of a solution to E g can be similarly estab- 
lished using Brouwer's fixed point theorem. We next study 
its uniqueness and the throughput optimality of a switching 
policy by resorting to an approximation given below, due to 
the complexity of E g . For the rest of this section, we will 
limit our discussion to the symmetric case where the channels 
have the same bandwidth and the system uses the same 
parameterization in all channels. We extend our discussion to 
more generic settings in the next section. 

Definition 5: A scheduling policy is unbiased if the station- 
ary channel occupancy distribution induced by such a policy 
is identical for every node, i.e., q^ = q^ for alH G Af and 
k G C. It is denoted by g u , and the space of unbiased policies 
is G u . 
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We can obtain an approximation to S § similarly as we did 
for E, using w w 



(fe) _ Q (k) Pi 



w+i 

2 



(a) 



A, 
~P 



1 



(J - 



t(i + e^)1} 



(b) 



and we have the following result. 

Theorem 3: Consider a system modeled by E § and the 
associated stability region A § . For all sufficiently large initial 
window sizes W, (i) the system of equations (E § , A) admits 
a unique solution, and (ii) g u is throughput optimal within 
the class Q u if = j± for all k. These are referred to 



as 



e qui- occupancy policies. 

Proof: We omit the proof on uniqueness, which is similar 
to the single-channel case; see Appendix [0 for the proof on 
throughput optimality. ■ 

The above results provide the following insights in addition 
to what we have observed in the single-channel case. Firstly, 
it's worth noting that E 8 reduces to E in the single-channel 
case by properly configuring related parameters, and E g thus 
constitutes a unified framework in describing the stability 
region of 802.11 DCF. 

Secondly, the uniqueness of the solution to (E § , A) is 
in fact true for even small windows. As an example, in 
Figure [3 we plot the numerical boundaries of stability regions 
for various window settings with equal channel occupancy. 
Compared to results in the single-channel case, convexity 
of the stability region is observed even with small backoff 
windows in the bi-channel case. Also, the numerical multi- 
equilibrium phenomenon disappears in this case. One way 
to explain this is by considering the discounting effect of 
channelization on the attempt rate. The attempt rate of each 
node in a channel is discounted by the occupancy probability 
in that channel. As discussed in the single-channel case, the 
attempt rate is roughly upper bounded by the reciprocal of 
the average backoff window size. Hence channelization has 
the effect of window expansion. The same explanation also 
applies to the observation that the stability region in a multi- 
channel system is nearly always convex. 
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Fig. 5. The stability region of bi-channel 802.11 DCF under the equi- 
occupancy policy. 
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Fig. 6. Throughput optimality of equi-occupancy distribution. 



Thirdly, given symmetric channelization, equal occupancy 
time is equivalent to equal packet assignment in each channel. 
The optimality of equi-occupancy policies therefore confirms 
the intuitive notion that load balancing (either in the number 
of active nodes or in the amount of date flow) optimizes 
the system performance in terms of expanding the stability 
region. In Figure we plot the analytical boundaries of 
stability regions corresponding to different unbiased policies in 
two scenarios. As can been seen, the equi-occupancy policy 
results in a stability region that is the superset of those of 
the other unbiased policies. It is also worth noting that as the 
backoff window increases, the gap between the superset region 
and other inferior regions decreases, as the reciprocal of the 
window size becomes the dominant factor in upper bounding 
the attempt rate. 

VI. Applicability and Implementation of Unbiased 
Policies in Both Symmetric and Asymmetric 
Systems 

In this section we discuss the applicability of the class of 
unbiased policies. We then present a number of practical im- 
plementations and their use in both symmetric and asymmetric 
systems. 

A. Unbiased policies 

We have so far restricted our policy space to unbiased 
policies that induce a node-independent channel occupancy or 
packet assignment distribution. Note that while nodes in the 
same system are typically programmed with the same protocol 
stack, the same protocol may not necessarily yield the same 
statistical behavior among different nodes. Nevertheless, there 
are a number of circumstances in which node-independent 
behaviors are induced, which justifies our focus on unbiased 
policies. Firstly, if the protocol explicitly prescribes packet 
allocation to each channel, the resulting packet assignment 
distributions are identical for all nodes. Secondly, if nodes have 
identical arrival processes, they then have unbiased behavior as 
well. Unbiasedness can also be observed in a saturated network 
(however, such a system is unstable). 

More generally, we note that when a node is active (i.e., its 
queue is non-empty and it is in the service process), from a 
mean-field point of view the channel conditions observed by 

(k) 

this node is fully characterized by p\ ' for each k (as a result 

(k) 

of the decoupling assumption), which is a function of r ■ ; for 
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all j ^ i. Therefore, the set of attempt rates {r^;Vz,Vfc} 
characterizes the contention condition in the system. If nodes 
are asymptotically symmetric, that is, liniAr^oo t„ 
for all i 7^ j and k, then we have 



k) i r (k) 



1, 



(k) 



lim ■ 

N^oo 2 



lim 

N^OQ 



i-n^(i-r/ fc) ) 



nw 1 



,(k) 



) 

Ar< fe) + (1-A) 



where A = riz^i j(l — r /^)- m m * s case we ma y consider 
the behavior induced by the underlying protocol on each node 
identical, and the corresponding policy unbiased. Note that the 
decoupling assumption is regarded as asymptotically true for 
a large number of nodes, so we may consider the asymptotic 
symmetry as an adjoint condition if we impose the decoupling 
approximation in modeling. 

B. Practical implementation of throughput optimal unbiased 
policies: symmetric channels 

We have shown that when channels are symmetric the op- 
timal switching policy within the class of unbiased policies is 
the equi-occupancy policy that balances load precisely. When 
channels are asymmetric, i.e., have different bandwidths, it is 
natural to expect that a load balancing policy yields throughput 
optimal performance, and to interpret a balanced load as 
having a packet assignment distribution proportional to the 
channel bandwidths. We will see that this interpretation is 
reasonable though not precise. 

We begin by commenting on how such policies may be 
realized in a symmetric system. 

We describe two very simple heuristics that implement an 
unbiased policy, and in particular, the equi-occupancy policy 
when channels are symmetric. The description is given in 
the bi-channel case for simplicity. The first is called SAS 
(switching after success), and the second SAC (switching after 
collision). In both schemes, a switching probability is assigned 
to each backoff stage. Under SAS (resp. SAC), a node switches 
to the other channel with probability ol{ upon a successful 
transmission (resp. collision) if it is in the ith backoff stage 
when this success (resp. collision) occurs. In addition, in SAC, 
after switching to the other channel, a node does not reset its 
backoff stage; instead, it continues the exponential backoff due 
to the last collision. Note that SAS can be used to implement 
any arbitrary packet assignment distribution (and thus load 
distribution), which is a useful feature when we proceed to the 
implementation under asymmetric channels. This is because 
with the assumption of nonpreemptiveness of the policy, i.e., 
A7, switching after each successful transmission is equivalent 
to assigning packets. 

These two schemes heuristically implement the equi- 
occupancy policy when channels are symmetric in the follow- 
ing sense. Consider the two-dimensional Markov chains for a 
bi-channel system in the form of Bianchi's model [7], where 
each state in one channel has a mirror state in the other. Since 
for both SAS and SAC, the corresponding Markov chain is 
irreducible with a finite number of states, using the argument 



of symmetry, the symmetric solution is the unique stationary 
distribution that reflects equi-occupancy. It should be noted 
however that neither of the above is a perfect solution and the 
key may be a proper combination of the two. The problem 
with SAS is that it can result in empty channels (the node that 
succeeded in the transmission happens to be the only node in 
that channel). When this happens nodes can tend to cluster in 
the non-empty channel for significant periods of time due to 
collision and backoff, while our mean field Markov analysis 
implicitly assumes no channels are empty for long. On the 
other hand, the problem with SAC (SAC rarely results in 
empty channels and avoids clustering in one channel) is that it 
interrupts the service process of a packet in a given channel, 
thus violating the nonpreemptive assumption about the policy. 

C. Practical implementation of throughput optimal unbiased 
policies: asymmetric channels 

We next proceed to asymmetric channels and examine how 
these two heuristics perform in this setting, and in doing 
so also empirically examine when the stability region is 
maximized. In particular, we focus on the performance of a 
policy when the majority of the nodes have similar arrival 
rates, and we examine the advantage of load balancing in 
improving stability. In our experiment, we fix 10 nodes with 
arrival rate 0.5Mbps that creates a mean-field background in 
a bi-channel network while inspecting the stability region of 
another two nodes, which is the projection of the aggregate 
stability region onto a plane of these two nodes' arrival rates. 
All nodes use the same policy in a single experiment. 
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Fig. 7. The projection of simulated stability region onto a plane of arrival 
rates of the two nodes under inspection. 

In Figure U\ we plot the empirical boundary of stability 
regions under different packet assignment distributions (im- 
plemented using SAS). As shown, policies with packet assign- 
ment ratio close to the bandwidth ratio indeed result in larger 
stability regions. However, while it seems safe to claim that 
properly balancing active time among channels according to 
their bandwidths improves the system performance, it remains 
unknown whether an exact match in load assignment is the 
optimal policy due to the nonlinearity of slot length in each 
channel w.r.t. active nodes. In addition, in practice we may 
not even know the effective bandwidth of each channel when 
channel conditions are imperfect. 

It is therefore highly desirable to have an adaptive mecha- 
nism that dynamically adjusts the load distribution in practical 
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implementation. Below we show that SAC to a large extent 
can achieve this goal, with the reason being that collision 
rate reflects the contention level and bandwidth information. 
Figure [7] also shows the empirical stability region obtained 
using SAC with switching probability at the zth backoff stage 
ai = m 9 w here m is the maximum backoff stage. SAC is 
clearly not optimal, but it maintains good performance under 
different bandwidth ratios. 
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Fig. 8. Histogram of node population in the slower channel: (a) ((b)) SAC 
(SAS) with on = 0.5; (c)((d)) SAC (SAS) with on = ^. 

We further highlight the adaptiveness of SAC in comparison 
to SAS. Assume that the active node population in each 
channel is the same and static, given then the same period of 
time, faster channels experience more transmission successes 
than slower ones. Therefore, if a SAS-like switching policy 
is adopted for a relatively congested network, nodes would 
cluster in the slower channels and the throughput performance 
degrades significantly. However, if the congestion is due to 
bandwidth asymmetry, then this is reflected in the collision rate 
of transmission, which in turns triggers channel reallocation 
under SAC. We illustrate this point using the following exper- 
iment. Consider a bi-channel system with strongly asymmetric 
channels, where the bandwidth of channel 1 (2) is 1Mbps 
(10Mbps). The system consists of 60 nodes each with an 
arrival rate 0.1Mbps, and this aggregate arrival rate (6Mbps) is 
slightly below the empirical saturation throughput under this 
setting. In the first test, we compare the resulting distribution 
of number of nodes in channel 1 between SAC and SAS with 
the switching probability = 0.5 for all stages in both 
channels, and we repeat the inspection with the switching 
probability oti = ^ at stage i in the second test; the duration 
of simulation is 180 seconds. The switching probability profile 
in the first test can be regarded as a blind configuration, while 
the second profile can be taken as an adaptive configuration 
that partially incorporates collision history into switching 
decisions. In Figure we plot the histograms of the number 
of nodes in channel 1, as well as the empirical throughput 



obtained. As can be seen, the blindly configured SAS drives 
nodes to cluster in the slower channel, while SAC avoids this 
problem. Interestingly, SAS has comparable performance as 
SAC if we adjust the switching probabilities as we did in the 
second test, which reflects the congestion level in the residing 
channel, and both distributions "match" the bandwidth ratio. It 
suggests that while SAS is not as adaptive as SAC, it remains a 
valid alternative implementation and could achieve comparable 
performance when configured appropriately, as did above. 

VII. Conclusion 

Using the characterization of the stability region of a 
multi-channel multi-user WLAN system, we investigated the 
throughput optimal channel switching schemes in such sys- 
tems. This work can be extended in the following directions: 
1) the effect of asymmetric channels on the characterization 
of stability region; 2) throughput optimal switching when 
considering the larger space of biased policies. 

Appendix A 
Proof of pi < pi 

We first define the following stochastic processes generated 
by the queueing process at node i. 

Ti i Q(t)/T i g(t) := the total length of real time periods up 
to time t that the queue at node i is 
non-empty /empty (or i is busy/idle); 
^,QW/^i qW := the t ^! number of slots up to time t 
that the queue at node i is 
non-empty/empty at the beginning 
of slots. 

These processes are well-defined on the same sample space 
Q. Assume that the queue is stable, then due to ergodicity pi 
and pi can be expressed respectively as 



Pi = lim — — = lim - 

*->°o t t^ooT itQ (u,t)+T ii Q(u,ty 



and 



pi = lim — 

t^oo N iiQ (u),t) + N if Q(<j,ty 

for all uo G ft. Let A^(t) be the total time fragmentation of 
busy periods in idle slots of node i up to time t, and let 
Si,Q(k) (S i Q(k)) be the length of the kth busy (resp. idle) 
slot. Quantities described above are illustrated in Figure |9] 
Then, we have 

T hQ (t) - Ai(t) = V # l0 (fc), 



k=i 



and 



Ni, Q (t) N i,QW 

t= £ S itQ (k)+ £ S.Q(k). 

k=l k=l 
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Fig. 9. Slotted time dynamics. 
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Therefore, 



> lim 



lim 

t— >-oo 



lim 

t— >-oo 



T itQ (t) - Ai(t) 



E£f (t) Si, Q {k)-\ 



E& (t) %^) 



^i,<?(*) 



iV jlQ (t) 



where we have suppressed the reference to a sample point 
uj in all involved processes for simplicity, or interpreted the 
equalities as with probability one. Let E[5^q] and E[S^ q] be 
the conditional average lengths of an arbitrary slot, given that 
the queue at node i is non-empty or empty at the beginning of 
slot, respectively. We claim that E[5^q] > E[5^g] (see the 
next Appendix for justification). Note also that Ni y Q(t) — >• oo 
and N { -q{1) —> oo as t —> oo due to the stability assumption. 
Consequently, following ergodicity, we obtain 

Ni, Q (t)E[Si, Q ] 



> lim 

- t ^oo N itQ (t)E[Si,Q] 



■ N- ■ 



fins* 



> lim 

t^oo 



N itQ (t) 



Ni,Q(t) + \Q(t) 



When the queue is unstable, we have pi = pi = 1. In either 
case, we have pi > pi. It remains to justify the claim made 
above, which appears in the next Appendix. 

Appendix B 

Computation of E [£{.}] and Related Quantities 

Given an event {•}, let P id i e] {-}, P S ucc- : {-} and P co u ]{ .} 
be the conditional probabilities that a slot is idle, that the 
transmission attempt in the slot is a success, and that the 
attempt is a collision, respectively. Notice that P co u-{.} = 



1 — Pidle;{-} 



P S ucc;{-}' Therefore, 



E[S { . } ] 



a-P 



idle;{-} 



■T„-R 



succ;{-} 



T r -P, 



coll;{-}- 



where a, T s and T c are the lengths of an empty system slot, 
a successful transmission, and a collision, respectively. Define 
then by Ti 5 Q the conditional probability that node i transmits in 
an arbitrary slot, given its queue is non-empty at the beginning 
of this slot, and hence we have r ? , 



Consequently, 



where 



Since Pa, 



Pidle;i,Q = (1 - T iiQ ) ]J(1 - Tj ) , 
Psucc;i,Q = ^2^i ]^[(1 - fj), 



"idie-,i,Q < Pidie;i,Q and a < min{T s ,T c }, we have 
E[S^q] > E[6' i g] and they are both finite. Explicit expres- 
sions for other variations of E [£{.}] can be derived in a similar 
way, and are thus omitted. 

When successive attempts are considered and the FACS 
decoupling approximation is adopted, we can adapt the above 
computation as follows. Denote by Z^dZe;{-}> L succ .^ and 
L co ii;{-} m e average lengths of the slot in the corresponding 
cases, and hence L icUe . { . } = a, L succ . { . } = T s and L coil . { . } = 
T c in the above computation. When successive attempts are 
taken into account, we have 



Lidle;{-} — 



J succ;{-} 



i=0 v J 1 w 



and 



L 



coll;{-} 



X 



CW 



{■} 



i \[\cw{ 

V cw { . } )\-± J 



CW 



{•} 



(cw { . } ) 
1 



(l + CW { . } )(l-£) 



1 ^S) 



where CW{.} is the conditional average backoff window size. 
These quantities are well-defined when W > 2 which is 
presumed in application. The first approximation of £ co ^ ; {.} 
is due to omitting the possibility of collisions involving three 
or more nodes, and the other one results from substituting 
CW{.} with the initial backoff window size W. Note that, if 
we neglect successive attempts, we have L succ .{.y = T s and 
L co ii-{-} = T c , which is also a natural consequence when W 
is sufficiently large in the above equations. 

Appendix C 
Approximation of fa 

Due to the analytical intractability of A^(t), we are inter- 
ested in proper approximations of pi that can lead to good 
estimate of A; a good estimate in the context of stability study 
means a tight underestimation. Recall that pi < pi and equality 
holds if and only if pi = 1 or pi = 0; therefore by replacing 
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pi by pi in E(c), solutions to the resulting system of equations 
form an underestimation of A but accurate when pi = 1 or 
for all i. Moreover, when < pi < 1, we have 



pi = lim 



T,, Q (£)-A,(t) 



(t)-A 4 (t) , T iiV (t)+Ai(t) 



< lim 



t^oo T,, Q (t) 



PiE[^] + (1 - Pi)E[5 i|Q ] 



< Pi 



where 



S£b(<) 



1 »,qW 



and defining 



we have pi < pi < pi. Hence, substituting pi with pi in E(c), 
we can obtain a tighter underestimation of A than with pi, 
thus trading off computational complexity for higher accuracy. 
Empirical results suggest that p is sufficiently close to p, and 
we use p as p throughout our computation. 

Appendix D 
Computation of Q c and Q p 

We first define the following processes generated by the 
channel switching policy at node i. 

(t) I 'N^ (t) := the total number of c-slots in channel k 
up to time t that node i is present at 
their beginning; 



(j) I S^- (j) := the length of the jth c-slot in channel k 
given the presence (resp. absence) of 
node i at its beginning; 
Ni(t) := the total number of n-slots at node i up to 
time t. 

When a channel switching is scheduled, the edges of c-slots 
of the two channels that a node switches between may not be 
aligned. Hence, there may be a period of unsynchronized time 
of the nodal backoff timer, as shown in Figure [TO) If omitting 
the unsynchronized time, we have 



E E s8o-)=* 

and then 

<>(*) 



N$(t) 



s ( + k hj) 



3=1 



E si -hi)> 



l^k j=l j=l 

lec 



channel 1 



channel 2 



! switch 



+ 



switch 



c-slot 



n-slot 



unsynchronized 
time 



Fig. 10. Illustration of channel switching and timer synchronization. 



or equivalently, 



,N«(t) (fe) 



Then can be expressed alternatively as 



q{ k ^ = lim 



= lim 



Ew^s^ l \t)/s^ l \t) 

lec y ' 1 



(k) , v / g^E[ftotgj] 



where 



JV<f>(t) 



We next define a few more generated processes at node i. 

P^ k \t) := the total number of packets transmitted by 
node i in channel k up to time t; 

B\ k \j) I Ii k \j) := m e number of busy slots in the 
service circle of the jth packet, 

where a service circle is the period in slots between the 
beginning of service processes of two successive packets. A 
busy slot refers to a slot in the service process and a idle slot 

(k) 

is a slot when the queue at the node is empty. Then, q\ J can 
also be expressed alternatively as 

(k) N™(t) 
q\ = lim 



lim 

t—too 



Aec lv +i 

ES )(t) (5f)(i) + 7f)(i)) 



EiecEjU {t \BP(j) + I?>(j)) 



lim 



where 



M (fe) 



E heC if°(*) 
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(*0, 



We then obtain 

J k ) - 



f_ i 



W(l-pW) 



zee #i ~ny 



(0 tf } 



r^(l-pP) 



where P„ 



Eg (t) Bf)q) 
p/ fc) (t) 



E[AT 



(fe) 



Appendix E 
Proof of Proposition 1 

Substituting E(b) in (a), we obtain 



2A; 



Ti = 



PQY + l) 
P(W + 1) 
P 



IF 1 

2 



2 ^ J 2 
A<((W-1)(7 + 2T) 



1 

t* + — —a + T 



P(W + 1) 



which can be rewritten as 



^((W- l)a + 2T) 
P(W + 1) 



Therefore, let y = £\ Tj , 7- = / (l + ^f?) and 
^1 + , and we have 



A,((W-l)a+2T) 
P(W+1) 



Then, E is equivalent to 



E : 



(a') 
(bO 



which admits only one solution, namely 



Appendix F 
Proof of Theorem 3 

Using E gt7 (a), we can rewrite E gt7 (b) as follows: 



p 



E W k) 

kec ^ 
f T( 1 + 



W- 1 



r (*o 



E-r)]} 

^Ef^E^f)^ 



where 0} = & = x ^ w -^ +2T 9 and 



MQ lk \pjJ * = 9 (fc) E^rj fc) = E^^-tgW] 2 with 



2Pi 
J — W+l 



> 



> for all j. Notice that (j>i{q^\ pj, j ^ i) is a 
convex function of given any fixed pj where j ^ i, and 
it is also an increasing function of p/s given any fixed q( k \ 
We then have 

* = 6£$>( g < fe >) + 0? 

kec 

where the equality holds when q^ = Therefore, when 
switching to the equi-occupancy policy from any arbitrary 
unbiased policy, the utilization factor of each node is always 
non-increasing. Hence, we conclude that the equi-occupancy 
scheduling policy is throughput optimal in Q u . 

Appendix G 
Miscellaneous 

Proof of the Markovian property of Qi(n) 

Let a n denote the length of the nth slot, and let {Pt}t>o 
be the transition semigroup of Qi(t). For any q i e N^i = 
1, 2, . . . , n, we have 

P(Qi(n + 1) = q n +\ I Qi(n) = q n , . . . , Q*(0) = q ) 

/ n+l n \ 

= PyQiC^Vj) = <7n+i I QiC^Vj) = q n ,---,Qi(0) =qo) 

/« / n+l n 
••• / =^n+l I QiC^tj) =g n ,..., 

Qi(0) = go,crj = tjj = 1, . . . ,n + 1^ x 

n 

X /a ri+ i,...,ai(Wl 5 ---^l I Q<(5Z^") = 4n,-..,Qi(0) = g ) 

J = l 

dt n +i ■ • • efti 



ri+i fe, Qn+l)fa n +i (^n+l) x 
n 

X fa nj ...,a 1 (tn, • • • , *1 | QiC^^j) = ?n, • • • , Q*(0) = ^o) 

dt n +i ■ • • d£i 

= ^ Pt n+1 (qn, Qn+l)fa n+1 (tn+l) dt n +\ X 
X J "' J /ffn,.,^^,---,*! I =9n, 



• • • ,Q*(0) = go) dt n • • - dti 
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J Pt n+ i(Qn,Qn+l)fcr n+ i(t 



n+l) 



where the third equality is due to the assumption that cr^'s are 
i.i.d. and independent of Qi(t), and / is the (joint) probability 
density function of cr^'s. On the other hand, for any g n , g n+ i G 
N, we obtain 

P(Qi(n + 1) = g n+ i | Qi(n) = g n ) 

✓ n+l n \ 

^ 3=1 3=1 ' 

/„ / n+l n 
••• / P[Qi(52t j ) = q n+1 | 

J V 3=1 3=1 

(7j =tjj = l,...,n + l^x 

n 

X fa n+1 ,...,a! (^n+l 7 • • • j ^1 | Qi( ^j) = <7n) ^n+1 ' * ' <^1 

3 = 1 

= j "j P t n +i(Qn,qn+l)f* n+1 (tn+l)x 

n 

X /o- n ,...,o-i(^n, ,^1 | Qi(^CTj) = <2n) ^n+1 " " * dt\ 

3 = 1 

= J Pt n+1 (qn,q n +i)fa n +i(tn+i) dt n +i X 

X J "' J /^n,...^i(W--^l | Qi(X^j) = 2n) 



dt n - -dt\ 



— J Pt n+1 (qn,qn+l)fa n+1 (tn+l) dt n +i. 

Consequently, Qi(n) is a Markov chain. 
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Total bandwidth 


11 Mbps 


Data packet length P 


1500 Bytes 


DIFS 


50 /js 


SIFS 


10 /xs 


ACK packet length (in time units) 


203 /is 


Header length (in time units) 


192 fjis 


Empty system slot time a 


20 /xs 


Propagation delay 5 


1 /xs 


Initial backoff window size W 


32 


Maximum backoff stage m 


5 


Data rate granularity A A 


100 Kbps 


Instability threshold constant 


1% 


Total simulated time Tf 


10 seconds 



TABLE I 

Specifications of the implementation of test bench. 
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